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Abstract 

We consider the problem of the statistics of the scattering matrix S of 
a chaotic cavity (quantum dot), which is coupled to the outside world by 
non-ideal leads containing N scattering channels. The Hamiltonian H of 
the quantum dot is assumed to be an M x M hermitian matrix with prob- 
ability distribution P{H) oc det[A2 + {H - e)2]-(/3M+2-/3)/2^ ^^^^^ ^ and 
e are arbitrary coefficients and /3 = 1,2,4 depending on the presence or 
absence of time-reversal and spin-rotation symmetry. We show that this 
"Lorentzian ensemble" agrees with microscopic theory for an ensemble of 
disordered metal particles in the limit M — > oo, and that for any M > it 
implies P{S) oc |det(l — ^t5'j|-{/3A'^+2-/3)^ where S is the ensemble average 
of S. This "Poisson kernel" generalizes Dyson's circular ensemble to the case 
S ^ and was previously obtained from a maximum entropy approach. The 
present work gives a microscopic justification for the case that the chaotic 
motion in the quantum dot is due to impurity scattering. 

PACS numbers: 05.45.+b, 72.10.Bg 
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Typeset using REVT^ 



I. INTRODUCTION 



Recent experiment sETq on conductance fluctuations and weak-localization effects in quan- 
tum dots have stimulated theoretical wor kin on phase-coherent conduction through cav- 
ities in which the classical electron motion can be regarded as chaotic. If the capacitance 
of the quantum dot is large enough, a description in terms of non-interacting electrons is 
appropriate (otherwise the Coulomb blockade becomes importantJlB). 

For an isolated chaotic cavity, it has been conjectured and confirmed by many examples 
that the statistics of the Hamiltonian H agrees with that of the Gaussian ensemble of 
random-matrix theory.^ If the chaotic behavior is caused by impurity scattering, the 
agreement has been established by microscopic theory: Both the Gaussian ensemble and 
the ensemble of Hamiltonians with randomly placed impurities are equivalent to a certain 
non-linear a-mo del.BS Transport properties can be computed by coupling M eigenstates of 
H to N scattering channels.0Tli Since N <^ M this construction introduces a great number 
of coupling parameters, whereas only a few independent parameters determine the statistics 
of the scattering matrix S of the system.0 

For transport properties at zero temperature and infinitesimal applied voltage, one only 
needs to know S at the Fermi energy Ep-, and an approach which starts directly from the 
ensemble of scattering matrices at a given energy is favorable. Following up on earlier work 
on chaotic scattering in billiardsjH§ two recent papersiil have studied the transport properties 
of a quantum dot under the assumption that S is distributed according to Dyson's circular 
ensemble.0'ii In Refs. ^ and ^ the coupling of the quantum dot to the external reservoirs was 
assumed to occur via ballistic point contacts (or "ideal leads" ). The extension to coupling via 
tunnel barriers (non-ideal leads) was considered in Ref. |lT[ In all cases complete agreement 
was obtained with results which were obtained from the Hamiltonian approach. iolli This 
agreement calls for a general demonstration of the equivalence of the scattering matrix 
and the Hamiltonian approach, for arbitrary coupling of the quantum dot to the external 
reservoirs. It is the purpose of this paper to provide such a demonstration. A proof of the 
equivalence of the Gaussian and circular ensembles has been published by Lewenkopf and 
Weidenmiiller,lli for the special case of ideal leads. The present proof applies to non-ideal 



leads as well, and corrects a subtle flaw in the proof of Ref. |T9| for the ideal case. 

The circular ensemble of scattering matrices is characterized by a probability distribution 
P{S) which is constant, that is to say, each unitary matrix S is equally probable. As 
a consequence, the ensemble average S is zero. This is appropriate for ideal leads. A 
generalization of the circular ensemble which allows for non-zero 5* (and can therefore be 
applied to non-ideal leads) has been derived by Mello, Pereyra, and Seligman, usmg a 
maximum entropy principle. The distribution function in this generalized circular ensemble 
is known in the mathematical literatureil as the Poisson kernel, 

P(5)oc|det(l-5t5)P^'""^-"\ (1.1) 

Here (3 G {1,2,4} is the symmetry index of the ensemble of scattering matrices: /3 = 1 or 
2 in the absence or presence of a time-reversal-symmetry breaking magnetic field; /3 = 4 



in zero magnetic field with strong spin-orbit scattering. (In Refs. ^ and |2^ only the case 



(3 = 1 was considered.) One verifies that P{S) = constant for 5" = 0. Eq. ( |1.1|) was first 
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recognized as a possible generalization of the circular ensemble by Krieger Ji for the special 
case that S is proportional to the unit matrix. 

In this paper we present a microscopic justification of the Poisson kernel, by deriving it 
from an ensemble of random Hamiltonians which is equivalent to an ensemble of disordered 
metal grains. For the Hamiltonian ensemble we can use the Gaussian ensemble, or any other 
ensemble to which it is equivalent in the limit M —* oo.il (The microscopic justification of 
the Gaussian ensemble only holds for M ^ oo.) For technical reasons, we use a Lorentzian 
distribution for the Hamiltonian ensemble, which in the limit M — »• cxo can be shown to be 
equivalent to the usual Gaussian distribution. The technical advantage of the Lorentzian 
ensemble over the Gaussian ensemble is that the equivalence to the Poisson kernel holds for 
arbitrary M > N, and does not require taking the limit M — oo. 

The outline of this paper is as follows: In Sec. ^ the usual Hamiltonian approach is 
summarized, following Ref. In Sec. |T|, the Lorentzian ensemble is introduced. The 
eigenvalue and eigenvector statistics of the Lorentzian ensemble are shown to agree with the 
Gaussian ensemble in the limit M — oo. In Sec. |I^ we then compute the entire distribution 
function P{S) of the scattering matrix from the Lorentzian ensemble of Hamiltonians, and 
show that it agrees with the Poisson kernel ( |1 . 1| ) for arbitrary M > N. In Sec. the Poisson 
kernel is shown to describe a quantum dot which is coupled to the reservoirs by means of 
tunnel barriers. We conclude in Sec. 0. 



II. HAMILTONIAN APPROACH 

The Hamiltonian approachEHllii starts with a formal division of the system into two 
parts, the leads and the cavity (see Fig. |T]a). The Hamiltonian of the total system is repre- 
sented in the following way: Let the set {|a)} represent a basis of scattering states in the 
lead at the Fermi energy Ep {a = 1, . . . , A^), with the number of propagating modes at 
Ep. The set of bound states in the cavity is denoted by (/x = 1, . . . , M). We assume 

M > N. The Hamiltonian Ti. is then given byEi 

n = Y: \a)EF{a\ + Y: \f^)H,Ai^\ + E {\f^)W,a{a\ + \a)W;M) ■ (2.1) 

The matrix elements if^j^ form a hermitian M x M matrix H, with real (/5 = 1), complex 
{(3 = 2), or real quaternion [jS = 4) elements. The coupling constants W^a form a real 
(complex, real quaternion) MxN matrix W. The NxN scattering matrix S{Ep) associated 
with this Hamiltonian is given by 

S{Ef) = 1 - 27riW\EF -H + mWW^Y^W. (2.2) 

For (3 = 1,2,4 the matrix S is respectively unitary symmetric, unitary, and unitary self-dual. 
Usually one assumes that H is distributed according to the Gaussian ensemble, 

P{H) = ^ exp (^-iPMX'hi H^) , (2.3) 

with V a normalization constant and A an arbitrary coefficient which determines the density 
of states at Ep- The coupling matrix W is fixed. Notice that P{H) is invariant under 
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transformations H — > UHU^ where U is orthogonal {/S = 1), unitary {(3 = 2), or symplectic 
(/3 = 4). This implies that P{S) is invariant under transformations W UW, so that it 
can only depend on the invariant W'^W . The ensemble-averaged scattering matrix 5* can be 
calculated analytically in the limit M — > oo, at fixed iV, Ep, and fixed mean level spacing 
A. The result is0 



It is possible to extend the Hamiltonian ( |2.1| ) to include a "background" scattering matrix 
5*0 which does not couple to the cavityS The matrix Sq is symmetric for (3 = 1 and can be 
decomposed as 5*0 = Oe^**0"^, where the matrix O is orthogonal and $ is real and diagonal. 
In the limit M — oo, the average scattering matrix S is now given byil 

S = (2.5) 

Lewenkopf and WeidenmiillerEl used this extended version of the theory to relate the 
Gaussian and circular ensembles, for /3 = 1 and 5 = 0. Their argument is based on the 
assumption that Eq. ( p. 51 ) can be inverted, to yield W'^W and S'o as a function of S. Then 
P{Si) = PsiS) is fully determined by 5* (and does not require separate knowledge of W^W 
and S'o). Under the transformation S USU^ (with U an arbitrary unitary matrix), S is 
mapped to USU'^, which implies 

PsiS)=PusuTiUSU^). (2.6) 

For 5" = one finds that P{S) is invariant under transformations S USU'^, so that P{S) 
must be constant (circular ensemble). There is, however, a weak spot in this argument: 
Equation ( |2.5|) can not be inverted for the crucial case S" = 0. It is only possible to determine 
W'^W, not S'o. This is a serious objection, since S'o is not invariant under the transformation 
S USU^, and one can not conclude that P{S) = constant for S = 0. We have not 



succeeded in repairing the proof of Ref. |T^ for S = 0, and instead present in the following 
sections a different proof (which moreover can be extended to non-zero S). 

A situation in which the cavity is coupled to n reservoirs by n leads, having Nj scat- 
tering channels (j = 1, . . . ,n) each, can be described in the framework presented above by 
combining the n leads formally into a single lead with = I]j=i Afj scattering channels. 
Scattering matrix elements between channels in the same lead correspond to reflection from 
the cavity, elements between channels in different leads correspond to transmission. In this 
notation, the Landauer formula for the conductance G of a cavity with two leads (Fig. |l|b) 
takes the form 

G=-^T. I^.P- (2.7) 

i=l j=Ni+l 
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III. LORENTZIAN ENSEMBLE 



For technical reasons we wish to replace the Gaussian distribution ( |2.3| ) of the Hamilto- 
nians by a Lorentzian distribution, 

where A and e are parameters describing the width and center of the distribution, and V 
is a normalization constant independent of A and e. The symmetry parameter (3 G {1,2,4} 
indicates whether the matrix elements of H are real [(3 = 1, Lorentzian Orthogonal Ensemble 
(AOE)], complex [/? = 2, Lorentzian Unitary Ensemble (AUE)], or real quaternion [f3 = 4, 
Lorentzian Symplectic Ensemble (ASE)]. (We abbreviate "Lorentzian" by a capital lambda, 
because the letter L is commonly used to denote the Laguerre ensemble.) 

The replacement of (p^ ) by ( ^TTP is allowed because the eigenvector and eigenvalue 
distributions of the Gaussian and the Lorentzian ensemble are equal on a fixed energy scale, 
in the limit M oo at a fixed mean level spacing A. The equivalence of the eigenvector 
distributions is obvious: The distribution of H depends solely on the eigenvalues for both 
the Lorentzian and the Gaussian ensemble, so that the eigenvector distribution is uniform 
for both ensembles. In order to prove the equivalence of the distribution of the eigenvalues 
El, E2, . . . , Em (energy levels), we compare the n-level cluster functions Tn{Ei, E2, . . . , En) 



for both ensembles. The general definition of the T„'s is given in Ref. |2^. The first two T„'s 
are defined by 

M 

Ti{E) = {J25{E-E,)) (3.2a) 

i=l 

M M M 

T2{Ei,E2) = { J2 ^{Ei - E,)6{E2 - E,)) - {Y.5{Ei - Ei)){Y.5{E2 - E,)). (3.2b) 

i,j=l i=l j=l 

The brackets (. . .) denote an average over the ensemble. The cluster functions in the Gaus- 
sian ensemble are known for arbitrary for the Lorentzian ensemble we compute them 
below. 

From Eq. (|3.1| ) one obtains the joint probability distribution function of the eigenvalues, 

Pm}) = iA^^(^^+2-^)/2 n \E^ - E/ n (A^ + (E, - ef))-'-'''''-'-'^^' . (3.3) 

i<j i 

We first consider the case A = 1, e = 0. We make the transformation 

S^\±^. (3.4) 

The eigenvalues e*'^-' of the unitary matrix S are related to the energy levels Ej by 

^i4>j ^ 1±1^ ^ = 2arctanE,-. (3.5) 
1 — lEj J J 
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The probability distribution of the eigenphases follows from Eqs. ( pl3|) and ( p.5|) , 

This is precisely the distribution of the eigenphases in the circular ensemble. The cluster 
functions in the circular ensemble are known.ciJ'c3 The n-level cluster functions in the 
Lorentzian ensemble are thus related to the n-level cluster functions in the circular 
ensemble by 

n 2 

T„^(Ei, . . . , = Tf (2 arctan^i, . . . , 2 arctan^) J] ^ ^- (3.7) 

i=i ^ + 

For n = 1 one finds the level density 

independent of (3. For n = 2 one finds the pair-correlation function 
rpA/j-, N 4 sin^(Marctani?i — Marctani?2) 

T2 (El, E2) = (1 + ^2)(i ^ ^2) sin2(arctanEi - arctan^a) " ^^"^^ 

Eq. ( ^.91 ) holds for [3 = 2. The expressions for /? = 1,4 are more complicated. 

The n-level cluster functions for arbitrary A and e can be found after a proper rescaling 
of the energies. Eq. (|3.71) generalizes to 



T„^(Ei, . . . , i?.) = Tf (2 arctan . . . , 2 arctan \{ (3.10) 

The large-M limit of the T„'s is defined as 

>;.(ei,---,en)= lim A"r„(eiA,...,e„A). (s.n) 

M^oo 

For both the Gaussian and the Lorentzian ensembles, the mean level spacing A at the 
center of the spectrum in the limit M ^ cxd is given by A = Avr/M. Therefore, the relevant 
limit M ^ 00 at fixed level spacing is given by M 00, A —>■ 00, A = Xtt/M fixed for 
both ensembles. Equation ( |3.10| ) allows us to relate the YnS in the Lorentzian and circular 
ensembles, 



>T(ei, • • • , a) = lim (27r/M)" Tf (2 arctan(7rei/M), . . . , 2 arctan(7ra/M)) 

M^oo 

= hm (27r/M)"Tf(27r6/M,...,2<n/M) 

M— >oo 



= F„^(ei,...,en). (3.12) 

It is known that the cluster functions Yj^ in the circular ensemble are equal to the clus- 
ter functions in the Gaussian ensemble.!! Equation (|3.12|) therefore shows that the 
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Lorentzian and the Gaussian ensembles have the same cluster functions in the large-M 
limit. 

The technical reason for working with the Lorentzian ensemble instead of with the Gaus- 
sian ensemble is that the Lorentzian ensemble has two properties which make it particularly 
easy to compute the distribution of the scattering matrix. The two properties are: 

Property 1: If if is distributed according to a Lorentzian ensemble with width A and 
center e, then is again distributed according to a Lorentzian ensemble, with width 
A = A/(A2 + e^) and center e = e/{\^ + e^). 

Property 2: If the M x M matrix H is distributed according to a Lorentzian ensemble, 
then every N x N submatrix of H obtained by omitting M — N rows and the corresponding 
columns is again distributed according to a Lorentzian ensemble, with the same width and 
center. 



The proofs of both properties are essentially contained in Ref. ^ In order to make this 
paper self-contained, we briefly give the proofs in the appendix. 



IV. SCATTERING MATRIX DISTRIBUTION FOR THE LORENTZIAN 

ENSEMBLE 

The general relation between the Hamiltonian H and the scattering matrix S is given 
by Eq. (|2.2|) . After some matrix manipulations, it can be written as 

^ = (l + inW^H - Ef)~^w) (l - zttW^H - Epy^Wy^ . (4.1) 

We can write the coupling matrix W as 

W = UQW, (4.2) 

where U is an M x M orthogonal {j3 = 1), unitary (/3 = 2), or symplectic {(3 = 4) matrix, 
14^ is an X matrix, and Q is an M x matrix with all elements zero except Qnn = 1, 
1 < n < A^. Substitution into Eq. ( |2.2| ) gives 

5 = (l + inW^HW) (l - irrW^HW) ~^ , (4.3) 

where we have defined H = Q^W{H - EpY^UQ. 

We assume that if is a member of the Lorentzian ensemble, with width A and center 0. 
Then the matrix H — Bp is also a member of the Lorentzian ensemble, with width A and 
center Ep. Property 1 implies that {H — Ef)~^ is distributed according to a Lorentzian 
ensemble with width A = A/(A^ + E^) and center e = Ep/{\^ + E'jp). Orthogonal (unitary, 
symplectic) invariance of the Lorentzian ensemble implies that U^{H — Ep)~^U has the 
same distribution as {H — Ep)~^ . Using property 2 we then find that H [being an A^ x A^ 
submatrix of U'^{H — Ep)~'^U] is distributed according to the same Lorentzian ensemble 
(width A and center e). 

We now compute the distribution of the scattering matrix, first for a special coupling, 
then for the general case. 
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A. Special coupling matrix 

First we will consider the special case that 

W = T^-^IH^^ (4.4) 
is proportional to the unit matrix. The relation (|4.3| ) between the S and H is then 



S - (4.5) 

Thus the eigenvalues Ej of H and e^'^^ of S are related via 

gi<^,. ^ = 2aTctanEj. (4.6) 

1 - lEj 

Since transformations H —>■ UHU^ (with arbitrary orthogonal, unitary, or symplectic 
N X N matrix U) leave P{H) invariant, P{S) is also invariant under S —>■ USU^ . So P{S) 
can only depend on the eigenvalues e**^^ of S. The distribution of the E's is [cf. Eq. (|3.3|)] 

pm) = 1a^(^-^+-^)/^ n 1^. - n (A^ + - ery''''^'-'''' . (4.7) 

j<fc j 



From Eqs. (|4.6|) and ( p7|) we obtain the probability distribution of the 0's, 

P({0,}) = 1 (i^) n - n |l - a^e^'i'' , (4.8a) 



j<k j 



1 - A - 2£ _ + - A - iE_ 
l + ~X + te ~ A2 + ^2 + A + 



^ = T = . J , , (4.8b) 



Eq. (O) implies that P{S) has the form of a Poisson kernel, 



det(l - SS^) 



{l3N+2~(})/2 



PiS) = " ciet(l - 5t^) , (4.9) 

the average scattering matrix 5* being given by 

Snm = CF5nm- (4-10) 



B. Arbitrary coupling matrix 

Now we turn to the case of arbitrary coupling matrix W . We denote the scattering 
matrix at coupling W by S", and denote the scattering matrix at the special coupling ( [4.4| ) 
by 5*0. The relation between S and 5*0 is 

S = r + t'Soil - r'So)-H ^ So = (t'y^S - r)(l - r^SyH^ (4.11) 
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where we abbreviated 



r = {1 - 7tW^W){1 + 7tW^W)-\ (4.12a) 

r' = ~W{1 - 7tW^W){1 + TTW^Wy^W-\ (4.12b) 

t = 2tt'^/'^W{1 + itW^W)-\ (4.12c) 

t' = 27c^^\l + 7rW^W)-^W\ (4.12d) 

The symmetry of the couphng matrix W is reflected in the symmetry of the 2N x 2N matrix 

Si=(lt], (4.13) 



which is unitary symmetric {j3 = 1), unitary (/? = 2) or unitary self-dual (/3 = 4). 

The probability distribution Pq of Sq is given by Eq. (|4.9|). The distribution P of S 
follows from 

P(S) = PoiSo)^, (4.14a) 

where the Jacobian dSoldS is the ratio of infinitesimal volume elements around Sq and S. 
This Jacobian is known,!!'!! 

dS l^|det(l-rt^)|2y ^^•■'^"^ 



After expressing 5*0 in terms of S by means of Eq. ( [4.11| ), we find that P{S) is given by the 



same Poisson kernel as Eq. (|4.9| ), but with a different S", 

l + ii{X + ie)Wm' 

In the limit M — »• cxd at fixed level spacing A = \ti/M ^ Eq. ( [4.15| ) simplifies to 

-_AM--K^W^W 

The extended version of the Hamiltonian approach which includes a background scat- 
tering matrix 5*0 can be mapped to the case without background scattering matrix by a 
transformation S ^ S' = USU^ {/S = 1), S ^ S' = USV (/3 = 2), or 5 ^ 5' = USU^ 
{j3 = 4), where U and V are unitary matrices.^ (U^ is the transposed of U, f/« is the 
dual of U.) The Poisson kernel is covariant under such transformations,^ i.e. it maps to a 
Poisson kernel with S' = USU^ (/? = 1), ^' = USV (/? = 2), or S' = USU^ (/? = 4). As 
a consequence, the distribution of S is given by the Poisson kernel for arbitrary coupling 
matrix W and background scattering matrix 5*0. This proves the general equivalence of the 
Poisson kernel and the Lorentzian ensemble of Hamiltonians. 
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V. IDEAL VERSUS NON-IDEAL LEADS 



The circular ensemble of scattering matrices is appropriate for a chaotic cavity which 
is coupled to the leads by means of ballistic point contacts ("ideal" leads). In this section 
we will demonstrate that the generalized circular ensemble described by the Poisson kernel 
is the appropriate ensemble for a chaotic cavity which is coupled to the leads by means of 
tunnel barriers ("non-ideal" leads). 

The system considered is shown schematically in Fig. 0. We assume that the segment of 
the lead between the tunnel barrier and the cavity is long enough, so that both the N x N 
scattering matrix Sq of the cavity and the 2N x 2N scattering matrix Si of the tunnel barrier 
are well-defined. The scattering matrix 5*0 has probability distribution Pq = constant of the 
circular ensemble, whereas the scattering matrix Si is kept fixed. 

We decompose Si in terms of N x N reflection and transmission matrices. 



Si= 7 J]- (5.1) 




The N X N scattering matrix S of the total system is related to 5*0 and 5*1 by 

S = ri+t[{l-Sor[)-'Soti. (5.2) 



This relation has the same form as Eq. (|4.11| ). We can therefore directly apply Eq. ( [4.14|) , 
which yields 

P{S) cx I det(l - rl5)|-(^^+2-/3)_ (5 3) 
Hence 5* is distributed according to a Poisson kernel, with S = ri. 



VI. CONCLUSION 

In conclusion we have established by explicit computation the equivalence for M > N 
of a generalized circular ensemble of scattering matrices (described by a Poisson kernel) 
and an ensemble of M x M Hamiltonians with a Lorentzian distribution. The Lorentzian 
and Gaussian distributions are equivalent in the large-M limit. Moreover, the Gaussian 
Hamiltonian ensemble and the microscopic theory of a metal particle with randomly placed 
impurities give rise to the same non-linear cr-mo Altogether, this provides a micro- 

scopic justification of the Poisson kernel in the case that the chaotic motion in the cavity 
is caused by impurity scattering. For the case of a ballistic chaotic cavity, a microscopic 
justification is still lacking. 

The equivalence of the Poisson kernel and an arbitrary Hamiltonian ensemble can be 
reformulated in terms of a central limit theorem: The distribution of a submatrix of H^^ of 
fixed size N tends to a Lorentzian distribution when M —>■ oo, independent of the details of 
the distribution of H. A central limit theorem of this kind for = 1 has previously been 
formulated and proved by Mello.ii 

This work was motivated by a series of lectures by P. A. Mello at the summerschool in 
Les Houches on "Mesoscopic Quantum Physics". Discussions with C. W. J. Beenakker, K. 
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APPENDIX: PROOF OF PROPERTIES 1 AND 2 OF SEC. Ill 



The two proofs given below are adapted from Ref. ^ The matrix H and its inverse 
have the same eigenvectors, but reciprocal eigenvalues. Therefore, property 1 of the 
Lorentzian ensemble is proved by showing that the distribution of the eigenvalues of is 



given by Eq. ( |3.3D , with the substitutions A — > A and e — > e. This is easily done. 



V . . 



e? 



-(/3Af+2-/3)/2 



1 ^MipM+2-p)/2 -Q \E,E,{Er^ - Ej^)f n 



V . . 

i<3 



d{E-' 



V . . 

t<3 



^(l3M+2-l3)/2 



1 ^M(/3M+2^/5)/2 -Q 1^-1 _ -Q ry2 ^ (^-1 _ .^,ymn2-P)/2 ^^^^ 



V 

In order to prove property 2, we may assume that after rescaling of H we have A = 1, 
e = 0. First consider M = N — 1. In this case, one can write 

H=(^^,iy (A2) 

where G is the N x N submatrix of H whose distribution we want to compute, F is a vector, 
with real (/? = 1), complex (/? = 2), or real quaternion elements (/3 = 4), and Z is a real 
number. For the successive integrations over Z and Y we need two auxiliary results. First, 
for real numbers a, b, c such that a > and 4ac > 6^, and for real m > 2 we have 

/oo ^ , — 7n+l/2 

dx {ax^ + bx + c)-"" = (ac - \b^) a"'~^7r^/^T{m - l/2)/r(m). (A3) 

-oo ^ ^ 

Second, if x is a (i-dimensional vector with real components, and if m > (d + l)/2, then 

POO 

dxi... dxd (1 + x^y^ = 7r''/2r(m - d/2)/r(m). (A4) 



oo 



Since det(l + H)^ is a quadratic function of Z, the integral over Z can now be carried 
out using Eq. (|A3|) . The result is: 



.2-j{-/3Af-2+/3)/2 ^ 



r . 7rl/2rfl(/?M+ 1-/?)) 

/ dZP(H) = ^ ^ det(l + 

J VT[^{PM + 2-/3)) 

(l + Ft(l + G^)-ly)-^^-^-^^ (A5) 
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Next, we integrate over Y. We may choose the basis for the F- vectors so that 1 + is 
diagonal, with diagonal elements 1 +Gf . After rescaling of the F- vectors to = Yi{l + Gf) 
one obtains an integral similar to Eq. ( [A4| ) , with d = j3{M — 1). The final result is 



{/3M-/3+l)/2 Y (l(3M +1-3)) 

Property 2 now follows by induction. Notice that Eq. ([A6|) allows us to determine the 
normalization constant V, 
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FIGURES 



FIG. 1. Schematic drawing of a disordered cavity (grey) attached to a lead. There are N 
scattering channels in the lead, which are coupled to M bound levels in the cavity. In (a) only one 
lead is drawn. A system with more leads (b) is described by combining them formally into one 
lead. 

FIG. 2. Schematic drawing of the chaotic cavity and the non-ideal lead containing a tunnel 
barrier. 
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Fig. 2 



